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Numerical instances are given of the speeding of the convergence of series by the Euler 
transformation. This is even applied advantageously to certain divergent series, and a 
rigorous justification is given. An example is given of a series for which use of the Euler 
transformation is not useful. Instances are given of several less widely known methods. 
Finally, the method of summation by transformation into a continued fraction is illustrated 
successfully in the case of certain divergent series. The possibility of applying two different 
methods in succession to a given series is exploited throughout the paper, in spite of the fact 
that this often requires summing a divergent series. 



A remarkably useful such transformation is based 
on a formula due to Euler: 

± a n x»=+- ± (A»a„) (V?-)"- (1) 

71 = 1 — X 71 = \1 — XJ 

We shall refer to this as Euler's transformation. 
A purely formal derivation of this is as follows. 
Recalling that 



a n = (l+A) n a 



we write 
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A rigorous derivation of eq 1 is given in [l], 3 pp. 
62-66. It is there stated that if the a's are all 
positive, then the right-hand side of eq 1 converges 
for every negative value of x for which the left-hand 
side converges. That this is not necessarily true 
for positive x is easily seen by taking x=l and 

1 



(n+l) 2 



Actually, it is clear that the right side of eq 1 is 
not likely to be an improvement over the left side 
unless x is negative and the a's are positive. In 
other words, the Euler Transformation might be 
expected to speed convergence only in the case of 
alternating series. However, for alternating seiies, 
the Euler Transformation is usually very helpful. 

The Euler Transformation is a special case of a 
more general transformation given by Markoff (see 
[9], pp. 178-194). However, we confine our atten- 
tion to the simpler Euler Transformation. 

1 The preparation of this paper was sponsored (in part) by the Office of Naval 
Research. 

2 Help in carrying out the computations was furnished by Nancy Mann and 
Shirley Marks under the direction of Gertrude Blanch. 

3 Figures in brackets indicate the literature references at the end of this paper. 



Applications of the Euler Transformation with 
x— 1 give 



log 2=J2 



(-1)- 



n 



« n2 n ' 



^=arctan 1 = 277 

4 #=6 2 n+l 



■l) n 



■ n=0 i- ' 



3-5- 



■(2»+l) 

(2) 



One recognizes the right side of eq 2 as the expan- 
sion of —log (1 — 1/2). 

An interesting variation of the Euler Transforma- 
tion is obtained as follows. 
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Applying this gives 
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Taking N=8 and x = l gives 



(N+l).»(N+m + l) 



log 2 = 1- 



2^3 



-i+i- 



1+1- 



1 



^2(9 2 9-10 



1 1-2 



1-2-3 



We get 



r+ 



4 9-10.11 ' 8 9-10-11. 12 



= 0.63452 381. 



Also (see table 1), the terms inside the curly brackets 
add up to 0.11724 674, and we get 

log 2 = 0.63452381+1 (.011724 674) = 0.69314718. 
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Table 1. 



1 
9 


=0.11111 111 


1 1 

2 9-10 


- 555556 


1 1-2 
4 9-10-11 


= 50 505 


etc. 


6 313 




971 




173 




35 




8 




2' 


0.11724674 



A more accurate value is log 2 = 0.69314 71806. 

One of the great advantages of the Euler trans- 
formation is that it can easily be used for numerical 
computations even in cases where one has only 
numerical values of the coefficients to work with; 
indeed, the computations arc very direct even in 
such cases, because only a simple differencing opera- 
tion is called for. 

A rather sensational use of the Euler transforma- 
tion is in connection with asymptotic series. We 
give an instance. 

We have for ?/>0 
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If the series within the square brackets were con- 
tinued indefinitely, it would diverge strongly. Thus 
it does not appear that the Euler transformation is 
applicable in the present case. However, let us pro- 
ceed formally to try to compute 

-be b Ei (—5) 

by use of the series alone. We get 



5^25 
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125 
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We have stopped just before the least term 24/625 
and if we add up the four terms shown the error is 
less than the next term. That is, we get 0.83200 
with an error less than 0.03840. Comparing with 
a more accurate value, 0.85211, we see that our 
error is indeed less than 0.03840. 

Formally, we have neglected the terms 



24 



120 



720 



625 3125 ' 15625 



and properly so, since they begin to diverge rapidly. 
Nevertheless, let us now try to sum them by use of 
the Euler transformation, taking x=— 1, 



a = 



a x ~- 



24 
= 625 

120 
= 3125 

720 
a2 ~15625 

etc. 

Wc readily form the differences of the a's numeri- 
cally (see the appendix for some 1 ricks for shortening 
the computation) and get 

a = 0.03840 00000 

Aa =0.00000 00000 

A 2 a =- 0.00768 00000 
A 3 a = 0.00307 20000 
A 4 a = 0.00645 12000 
A 5 a =0.00761 85600 
A 6 ao = 0.0140() 97600 
A 7 a = 0.02602 59840 



A s a = 0. 05613 40416 



A 9 a =0. 13145 60410 
A 10 a =0.33766 21486 
A%o=0.93823 63791 

A 12 a -2.80697 67610 
A 13 a =8.98851 15364 
A l4 a =30.66826 95734 
A 15 a ft = 111. 03898 71074 
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So the Euler transformation, if applied formally, 
gives the series 

10- 10 {1920 00000 + 96 00000-19 20000 + 20 16000 

-11 90400 + 10 99200-10 16640+10 96368 

-12 83750 + 16 48741-22 90616 + 34 26485 

-54 86152 + 93 59213-169 43205 + - • •} 

This is actually a divergent series, and the last 
terms listed are already increasing. Nevertheless we 
treat it like an asymptotic series, and add up all 
terms preceding the one that is numerically least. 
This gives 0.02016 04800. Adding this to the 
0.83200 that we had already, we get 0.85216 04800, 
which compares favorably with the more accurate 
value 0.85211. 

We neglected the terms 

-10- 10 {10 16640-10 96368 + 12 83750-16 48741 

+ 22 90616-34 26485 + 54 86152-93 59213 

+ 169 43205-- • •} 

We now apply the Euler transformation to these 
terms, taking x= — l, and a = 10 16640, a 1 = 1096368, 
• • •. This gives 

-10" 10 {5 08320-19932 + 13457-4372 + 916-1383 

+ 50-570 + 6 ■ • •} 

This adds up to -0.00004 96486 if we stop before 
the least term. Adding this to our previous value 
gives 0.85211 08314, which compares very favorably 
with the more accurate value, 0.85211 08814. 

The agreement with the more accurate value is 
experimental evidence that use of the Euler transfor- 
mation is justifiable even in the present case where 
the series diverges. Actually, we can justify use of 
the Euler transformation rigorously for the present 
series. Looking at eq 4, we see that the remainder 
can be written 
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which can be expanded as 
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If we put N=3 and y = 5, we easily verify that 

We then see that all differences of a Q of even order 
are positive, and that if we stop with a term involving 
a difference of odd order, the error involved is less 
than the next term (which involves a positive differ- 
ence of even order). 

To justify a second application of the Euler trans- 
formation, we modify the remainder term given 
above as follows: 
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We can apply the Euler transformation to other 
asymtotic series with equal success. For instance 
w T e have 
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If we now write 
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and expand 
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we can justify applying 
the asymptotic series 



Eider transformation to 



1 



1-3 



2w 2 ' (2w 2 ) 2 



We applied the Euler transformation to this scries 
with w=Q to compute a 20 decimal value of 



12e 36 f°° e-x 2 dx, 



which was then checked by a different means of 
computation (see [3], p. 86). 

We can use the Euler transformation to derive 
Airey's "Converging Factor" (see [2]). In those 
cases (as above) in which we can justify using the 
Euler transformation on an asymptotic series, we 
are then afforded a means of justifying the use of 
Airey 's "Converging Factor." 

To see how this is done, let us write. 
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Airey refers to the series in the curly brackets as 
the "Converging Factor" and gives a purely formal 
device for estimating a value of it (Airey 's estimate 
is justified in [12]). We will now show how to derive 
Airey 's formulas by use of the Euler transformation, 
which we have already shown to be justified for the 
present series. We write 

Then the "Converging Factor" takes the form 

We now apply the Euler transformation with 
as=-/S, 
a„—l 



-O-'X^) 



etc. 



We find 
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etc. 
So use of the Euler transformation gives us 
1 ( . , _J_ h ( <3 V (\ h(\-h )\ 

+(TT-J(^ L '- M ^^) + -i 

If we rearrange this in powers of 1/p (such a rear- 
rangement needs justification) , we find that it is 

identical with the formula on p. 526 of [2], 

If we choose v=N, y=N, so that h = 0, = 1, the 
"Converging Factor" reduces to 



2\ l+ l£i(iN) H ) 



(5) 



with Ci=l f C 2 = -1, Gi= -1, C 4 =13, etc. The 
values of the first 22 CVare given on p. 529 of [2]. 
Using the first 17 of these with N=^5 we get the 
following series for the "Converging Factor": 

i 10" 15 {1 00000 00000 00000 + 5000 00000 00000 
-250 00000 00000- 12 50000 00000 + 8 1250000000 
-1 46875 00000-11406 25000+19117 18750 
-6566 79688-303 73047+1686 04590-927 48286 
+ 26 80585 + 364 18569-295 53896+41 19945 
+ 154 50012-177 02724+ . . . }. 

This series is probably divergent if carried far 
enough, but we sum it as if it were an asymptotic 
series and stop just before the least term shown. 
Thisgives the value 0.52372 086. Airey gets 0.52372 087 
(sec [2] p. 529), but the difference of 1 unit in the 
last place is likely due to differences in rounding pro- 
cedure. This should be multiplied by 0.0384 and 
added to 0.832, which gives 0.85211 08810 as com- 
pared with the more accurate value 0.85211 08814. 

In [2], Airey applies the "Converging Factor" to 
various asymptotic series without rigorous justifiea- 
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tion, but with great numerical success. By deriving 
the "Converging Factor" by means of the Euler trans- 
formation, one could justify its use in some of the 
cases given in [2], though there would still remain 
the problem of justifying various rearrangements 
made by Airey of the series in the "Converging Fac- 
tor." Some of these can be justified by the pro- 
cedure in [12]. 

Although we have seen that the Euler transforma- 
tion can be applied in some rather remarkable situa- 
tions, nevertheless there are cases where it cannot 
be applied. In particular, for non alternating series 
(for instance, when x and the a's are all positive in 
(1)), it is of almost no use. However even for alter- 
nating series it occasionally fails, and we will now 
give an alternating series for which it is of no value 
whatever. 

Let us expand 



I 



dx 

1+x 2 



by the Euler-Maclaurin sum formula (see [4], p. 128) 
using an interval of length unity. We first need the 
values of the derivatives of 



1+x 2 



of odd order evaluated at and 1. At 0, the deriva- 
tives of odd order are zero. To get the values at 1 
we write 
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Putting x = l and recalling 
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Substituting in the formula 



l F(x)dx=±F(p)+±F(l) 
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(see [4], p. 128, with o=l, r=l) we get 
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73 54368 5 65248 ' 98304 
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This is certainly an alternating series. However, 
the partial sums do not lie alternately above and 
below the true value. In table 2 we have given the 
terms of the series with the partial sums to the right. 
If we try to improve the value bv use of the Euler 
transformation, we get 10- 5 (79U7+{12 — 10+61 — 
758 + - • •}), which is certainly no better. 



Table 2. 


Series 


Partial sums 


+0. 75000 


0. 75000 


+4167 


. 79167 


-50 


.79117 


+24 


. 79141 


-65 


. 79076 


+596 


. 79672 


-0. 13743 


. 65929 


+6. 69282 


7.35211 


True vali 


ie 0.7854 



As indicated, the Euler transformation is valueless 
for speeding convergence of a nonalternating series. 
Occasionally one can use a trick to transform a non- 
alternating series to an alternating series. Thus the 
formula 

n ={ n 71=1 n 
enables us to get a value for 

CO 1 

t^i n 2 

by (indirect) use of the Euler transformation. How- 
ever, in general such tricks are not available. 

Fortunately the Euler-Maclaurin sum formula is of 
great value for summing many nonalternating series. 
For instance, by use of the Euler-MacLaurin sum 
formula (see [4], p. 128) we prove that for ra> 1, s>l, 
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m 1 " l 1 
?? = i rl> n-1 to 
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The scries on the right is divergent, but it is an 
asymptotic series, and by taking m large compared to 
s, considerable accuracy can be obtained, which can 
be further increased by applying the Euler trans- 
formation to the asymptotic series. Thus, consider 
the case where 5=2. We have the well known result 



£{ n 2 6 



= 1.64493 40668. 



Even with m as small as 2 we get 
1,1 1 



1+V- 
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691 



3617 



43867 



223 64160 ' 1 9G608 668 46720 ' 4183 81824 
1 74611 8 54513 



6920 60160 ' 11576 27904 



This is divergent, but if we sum all terms before 
the least, which is 

, 691 



223 64160 



we get 1.64495 as compared to the more accurate 
value 1.64493. If we now apply the Euler trans- 
formation to the terms neglected, we get the series 

-10- 10 {1 54488-11765 + 17249-11523 
+ 14389-20493+ . . . }. 

Doubtless this diverges, but if we sum all terms 
before the least, we get 1.64493 36 as compared with 
the more accurate value 1.64493 41. Another use of 
the Euler transformation on the terms neglected 
above gives further improvement. 

In [5], Bickley and Miller give a method for deal- 
ing with series Xa n in which 



a n -i 



n n 2 n 3 
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In those cases in which a n is a differentiable func- 
tion of n, one can treat such series by the Euler- 
Maclaurin sum formula, and usually one gets the 
same approximation as would be given by the 
Bickley-Miller method . 

In the Bickley-Miller method, one writes 



n=0 



S n 



=5j a n- 
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Then one tries to determine constants a_,, a , a u 
... so that 



S-S„ 



m in 2 
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In general, this is impossible, since the series 



i i a i i 
m 



will be divergent for each in. However, usually this 
series is an asymptotic series, so that for large in- 
great accuracy is attainable. 

To determine the a's, we proceed formally as 
follows. By eq 7, 



S — S m -i = a m -A a-\(m — l) + a 
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Subtracting eq 7 from this gives 
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However, by eq 6, 
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Substituting this arid factoring out a m _i and dividing 
by in gives 



1 
in m 
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If we now expand 

1 1 r r(r+l) 
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and equate powers of — ; we can determine the a's. 
m 

Indeed we have 

a_ 1 (/l 1 -l)=l, 
(ao+1) A^a-xAz, 

a l (A l +l) = (a +l)A 2 +a-iA 3 , 
etc. 

(Note error in sign in formula (12c) on p. 757 of [5]). 
If we apply the Bickley-Miller method to the series 

ao 1 

si 
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we get the same expansion that is given by the Euler- 
Maclaurin sum formula (see above). Bickley and 
Miller do not apply their method to any series of 
interest that cannot be handled by the Euler-Mac- 
laurin sum formula, but they indicate how such series 
might arise. Bickley and Miller do not give a rigor- 
ous justification of their method, but Szasz (in [10]) 
has discussed the conditions under which one can use 
at least the term a.„ x m in eq 7 for a first approxima- 
tion to S. It would be very helpful to have a general 
treatment of the conditions under which the Bickley- 
Miller method can be used. 

In [7] and [8], Stieltjes gives some interesting ways 
for estimating the error in certain well known 
asymptotic series. However, his methods are rather 
special, and can hardly be applied to series at random. 

We turn now to what is perhaps the most widely 
applicable method that we know. This is the 
method of transforming the series into a continued 
fraction. 

Let us return to the example given above in which 
we got an approximation for 

Jfci n 2 

by the Euler-Maclaurin sum formula. We had a 
divergent series, of which we neglected all terms 
beginning with the least. Let us now try to find a 
sum for these neglected terms by transforming them 
into a continued fraction. To write them as a 
continued fraction, we first consider the series 
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Transforming this into a continued fraction by the 
procedure given on pp. 196 to 202 of [6] we get 

Oq 0]_ 02 O^ 

Z+Z+Z + Z + • • • 

with the coefficients 

a =0.0000 30897 65053 
a 1 = 1.1523 15485 
a 2 =0.3674 32414 
a 3 = 1.7288 91714 
a 4 =0.8615 52894 
a 5 =2.4273 31997. 



(Our a's after a are taken to be the negatives of those 
defined in [6]. Due to a gradual loss of significant 
figures, the final digits of our a's are increasingly 
inaccurate.) 

We then put 2=1, and use the resulting continued 
fraction to compute an estimate for our remainder 
(which, by subtraction from the known sum of the 
series should come out to be —0.00001 54985). The 
successive convergents are: 

-0.00003 08977 

-0.00001 43555 

-0.00001 67677 

-0.00001 53294 

-0.00001 57006 

-0.00001 54620. 

We note the characteristic continued fraction be- 
havior of the convergents in that the}^ are successively 
above and below the true value, and mono ton ely 
approach the true value from above and below. 
Already the sixth convergent gives a value correct 
within less than 4 in the eighth decimal place, and it 
seems very plausible to conjecture that later converg- 
ents will continue to come ever closer, in spite of the 
fact that it is a divergent series that we are trans- 
forming into a continued fraction. This is a not 
uncommon phenomenon with continued fractions; in 
fact this is the basis of the Stieltjes method for sum- 
ming asymptotic series. In [6] on pp. 349 to 361, p. 
365, p. 367, and pp. 372 to 373 are given many 
examples of convergent continued fractions, which 
are the formal equivalents of divergent asymptotic 
series, including some that resemble very closely the 
series with which we are dealing. 

Note that we used only odd powers of z in our series 

691 1 

223 64160 z ' ' 

This has the effect of greatly reducing the labor of 
finding the corresponding continued fraction (in the 
notation of pp. 196 to 202 of [6], all the b's will be 
zero) . 

For an additional instance of the use of continued 
fractions, let us recall the series for the "Con verging 
Factor" which we had earlier. We used all terms 
preceding the least, leaving several terms unused. 
Let us now try to attach a value to these unused 
terms by transforming them to a continued fraction. 
To avoid difficulties with the continued fraction 
expansion, we start with the term after the least. 
We formally transform the series 



10" 



364 18569 295 53896 , 41 19945 



154 50012 177 02724 
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into a continued fraction 



getting 





a ai 

2 + 2 +"" 


flo= 


0.00000 00182 09284 


fli= 


0.81150 624 


a 2 =- 


-0.67210 178 


as = 


0.80673 455 


a 4 = 


0.02419 193. 



Because one of the a's is negative, we no longer get 
the familiar performance of the convergent^. Never- 
theless, they are apparently converging to the true 
value 0.00000 00078 69 since they are 

0.00000 00182 09 

0.00000 00100 52 

0.00000 00052 40 

0.00000 00079 44 

0.00000 00079 16 

respectively. Using the last, we get a value of 
0.85211 08814 25 for —5e?Ei(—5) as compared with 
the more accurate value 0.85211 08814 237. 

One other remarkable use of continued fractions 
is given in [11], where the asymptotic series 



i 



1 



1 



ey2(ly 2w + 4w 3 + 8w 5+ 16w 7 



1-3 



1-3-5 , 1.3.5.7 



32w 9 



is converted into the continued fraction 

_i_ JL A JL _§_ 
2w — 2w — 2w—2w—2w— 

This continued fraction was used to compute por- 
tions of a table (see Section 26 of [1 1]) in spite of the 
fact that the continued fraction is likely divergent 
due to the fact that occasional convergents will have 
extremely small denominators. However, as shown 
by the estimates in Section 8 of [11], whenever we 
have a convergent with a large denominator, this 
convergent is a close approximation to the value of 
the function. 

In closing, we might remark that some of the 
transformations studied in theories of summability 
are occasionally useful to speed convergence. Thus 
the use of Riesz sums, 



&n~' 



2_j X m O „ 

m = l 

n 



where the \ m are increasing positive numbers, is 
sometimes helpful in estimating the sum of mi alter- 
nating series of which S m are the partial sums. 

Appendix. Methods of shortening the 
computation of A m a . 

In an earlier draft of this appendix, we set forth a 
scheme for shortening the computation of A m a . 
Upon reading this earlier draft, Prof. J. C. P. Miller 
proposed the following very ingenious scheme for 
computing A m a . 

In the text, we put 



a = 



d\ = 



24 
625' 



120 



a 2 ~- 



3125 

720 

I5()25' 



etc. 



and wished to compute A m a . To do this more 
quickly, we write 

n _ _ (N- i)l 



where N=b and 



«o=l' 



«!= l> 



o»=l + 



N' 



«H 1+ K>( 1+ n)> 
*-( 1+ y)( 1+ ^)( 1+ ? 



etc. 



Then 



(N—l)\ 



'(-N) N ~ l 
So we seek a formula for A m a n 
Actual subtraction gives 



Aa n =jja n . 
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Then we prove by induction on m that 



Putting n=0 gives 



AT 



A m ^a n = '^^A m a n +^A m - l a n . 



m t 



A^^o^lA^o+A^^o). 



(A) 



Then 



m , 



A w+1 a =-^(A-ao+A w - 1 a ). 



Recalling that it is actually 



A m a, 



that we seek, we write the simple recursion 



*+i 



a n m 



( A m a A m - l a \ 

I ^ om + l "I O w / 



2^+2 4AT 
Using this, we readily compute numerical values of 

A m a 



In our earlier scheme, we had expressed A m a as 
a polynomial in (1/iV), and listed the following re- 
sults (which can easily be verified by use of eq A) : 

Aa = 



A 2 a = 



A 3 «o = ^> 

_ 6 _i_ 3 

Aa °~N z ^~N 2 ' 
- _24 , 20 

120 . 130 ,15 

720 . 924 , 210 

5040 . 7308 ,2380 105 



40320 . 64224 , 26432 , 2520 



3 62880 , 6 23376 , 3 03660 , 44100 , 945 



A 11 « = 



36 28800 , 66 36960 , 36 78840 . 7 05320 , 34650 



JSfio 



h 



N* 



+ 



N* 



r 



N 7 



+ - 



N* ' 



^ 399 16800 769 98240 , 473 24376 110 98780 

A -Q; Arii \rm atq 



flu 



iVio 



iV9 



iV 8 



8 66250 10395 



4790 01600 , 9675 24480 . 6475 36032 

Al3 <*0= A712 4- ATT1 + " 



AT12 



N u 



N 10 



1773 31440 188 58840 5 40540 

A.79 ' A7"8 



iV 9 



N* 



N 7 



62270 20800 , 1 30967 36640 , 94189 45536 



AT13 



#12 



AT11 



29205 25608 3894 49060 182 88270 1 35135 

Ann 7vra ~> at-r 



ATio 



N» 



AT8 



N 7 



8 71782 91200 , 19 00603 35360 , 14 54105 80224 

A lo «0 = Am + " 



JSfu 



AT13 



ATi2 



4 9 9528 6296 79349 27000 , 5200 59540 9-1 59450 . 

ATU A710 ATQ 



N n 
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